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ABSTRACT 


In  this  paper  we  consider  an  approximation  scheme  for 
an  optimal  control  problem  described  by  a  hyperbolic  partial- 
functional  differential  equation  used  to  model  the  elastic  motion 
of  a  viscoelastic  body  of  Boltzmann  type.  The  method  is  based 
on  combined  finite  element/averaging  approximations.  We  present 
theoretical  and  numerical  results  for  a  problem  with  quadratic 
cost  functional. 


Introduction 

We  discuss  the  development  of  numerical  schemes  for  the  modeling  and 
control  of  longitudinal  vibrations  in  a  rod  with  Boltzmann  -  type  viscoelastic 
damping.  It  has  been  thought  for  some  time  that  a  high  fidelity  model  with  an 
accurate  description  of  the  damping  mechanism  is  crucial  for  applications  in¬ 
volving  control  and  stabilization  of  large  flexible  structures.  The  model  which 
we  describe  below  may  be  considered  as  a  prototype  for  the  investigation  of 
such  applications.  In  particular,  we  consider  Boltzmann-type  viscoelastic 
damping  in  our  model  (other  types  of  damping  which  have  been  studied  include 
viscous,  Kelvin-Voigt,  structural,  hysteretic,  etc.). 

Our  ideas  apply  and  extend  to  various  structures,  but  for  simplicity  of 


Si 

*  r 


exposition  we  consider  a  model  for  the  longitudinal  motion  of  a  uniform  bar 
with  fixed  ends.  This  leads  to  the  following  equation  (see  [5]): 

Mil  „  a2  a  f  a  f°  a 

P  ___u(t,x)  =  _.«_u(t,x)+  g(s)  — u(t+s,x)ds  ■+  f  (t,x) 
dr  ox  dx  dx 


where  p  is  the  density,  and  f  is  the  (applied)  body  force.  Here,  u(t,x)  is  the 


longitudinal  displacement  at  position  x  along  the  rod  at  time  t.  The  constant  a  is  a 
stiffness  parameter  and  the  function  g(s)  (described  in  more  detail  below)  may  be 


considered  as  the  damping  parameter.  Note  that  without  the  integral  term  this 
would  be  a  purely  elastic  model  (the  wave  equation  with  no  damping).  The 
integral  term  (referred  to  as  the  memory  or  history  term)  arises  from  the  fact  that 
in  the  underlying  constitutive  equation  for  the  Boltzmann  model,  the  stress  is 
assumed  to  be  a  function  of  the  strain  and  the  strain  history. 


l7>'t 


We  make  the  following  reasonable  physical  assumptions  on  the  function 


g:  (-“,0]-JR  and  constant  a: 


(1)  g  €  H^-r.O),  a  >  0. 


(2)  There  exists  a  continuous  function  g0  :  (-*,0]  -<-“,0] 

and 


constants  ti  >  0,  €n  >  0  such  that 


(a)  g(s)  =  g0(s)  =  0,  i  t  r 


(b)  g(s)  S  g0(s)  <  0,  -r  <  »  S  0 

(c)  €  i  a  +  |°r  g(s)ds  >  £0 


(d)  —  g(s)  <  Mg0(s)  a.e.  on  [-r,0]. 

ds 


The  existence  of  £n  follows  from  general  properties  of  elastic  moduli. 


Condition  (d)  is  a  "decaying  memory"  assumption.  For  a  further  dis¬ 


cussion  of  the  physical  basis  for  these  assumptions,  see  [15],  [16],  and 


[24].  The  fixed  end  boundary  conditions  are  given  by 


u(t,0)  =  0  =  u(t,i). 


We  consider  initial  data  of  the  form 


u(0,x)  =  d(x),  —  u(0,x)  =  v(x),  0  <  x  <  1 

ot 


— u(s,x)  =  h(s,x),  -r  <  s  <  0,  0  <  x  <  1. 

ox 


This  represents  initial  displacement,  velocity  and  past  history. 
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We  are  interested  in  the  problem  of  constructing  a  sequence  of  finite 
dimensional  "approximate"  models  which  can  be  used  for  control  design.  Our 
ideas  are  developed  in  the  context  of  standard  results  from  linear  semigroup 
theory,  such  as  the  Hille-Yosida  theorem  and  the  Trotter-Kato  approximation 
theorem.  In  addition,  we  appeal  to  many  results  concerning  the  abstract  linear 
quadratic  cost  optimal  control  problem  ([12],  [13],  [7]).  Again  we  note  that  our  ideas 
can  be  easily  modified  to  include  other  boundary  conditions,  as  well  as  easily 
extended  to  include  structures  such  as  Euler-Bernoulli  beams  with  Boltzmann 
damping. 


In  section  2,  we  develop  a  state  space  formulation  of  the  class  of  PFDE’s 
which  we  will  consider.  A  well-posedness  result  is  given  in  this  context.  In 
section  3,  we  develop  an  approximation  scheme,  and  convergence  results  are  given. 
In  section  4,  we  present  some  examples,  including  an  application  to  a  quadratic  cost 
optimal  control  problem. 


The  notation  in  this  paper  is  standard.  The  symbols  <.,.>x  and  H * H x  stand 
for  the  inner  product  and  norm,  respectively,  on  the  Hilbert  space  X.  Often 
the  underlying  space  X  is  not  specified  but  will  be  understood  from  the  context. 
The  symbol  =  indicates  that  the  expression  on  the  left  is  defined  by  the 
expression  on  the  right.  Also,  E(X,Y)  denotes  the  space  of  bounded  linear 
operators  from  X  to  Y. 
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2.  A  State  Space  Model 


We  now  proceed  with  the  formulation  of  the  system  (1.1 )-( 1 .3)  as  an 
equivalent  abstract  Cauchy  problem  in  an  appropriate  Hilbert  space.  In  several 
recent  and  relevant  investigations  ([10],  [11],  [22],  [23]),  the  underlying  PFDE  is 
treated  as  an  abstract  FDE.  In  a  manner  analogous  to  ordinary  FDE’s,  a 
solution  semigroup  is  defined,  its  infinitesimal  generator  is  characterized,  and  a 
Cauchy  problem  is  formulated  on  a  "product"  space.  An  important  problem  in 
this  approach,  however,  is  the  prescription  of  appropriate  initial  data  so  that  the 
abstract  FDE  is  well-posed.  Kunish  and  Schappacher  have  shown  ([18])  that  a 
"natural"  choice  for  a  product  space  generally  leads  to  a  Cauchy  problem  which 
is  not  well-posed.  The  correct  choice  of  state  space  often  involves  the  use  of 
suitable  interpolation  spaces  (sec  [9],  [10]). 

We  proceed  in  ?  different  manner  by  first  defining  the  state  space  (a 
"product"  space)  Z  and  state  operator  A  .  Wcll-posedness  follows  when  it  is 
shown  that  A  is  the  infinitesimal  generator  of  a  C0-semigroup  on  Z  .  To 
proceed,  it  will  be  convenient  to  introduce  the  space  L°  =  L°(0,1)  defined  by 


j<p  e  L2(0,1)  :  |  <Kx)dx  =  0 


Lj  =  Lj(0,l)  =  <p  €  L2(0 


To  explain  the  use  of  the  space  L°,  we  remark  that  (as  will  be  described  below) 
we  choose  to  formulate  the  second  order  (in  time)  differential  equation  (1.1)  as  a 
first  order  system  using  strain  ux(t,x)  and  velocity  ut(t,x)  as  states  rather  than  dis¬ 
placement  u(t,x)  and  velocity.  (This  is  done  because  we  use  strain  and  velocity 
feedback  in  the  optimal  control  problem  to  be  considered  below  in  Section  3).  The 
space  L°  corresponds  to  that  part  of  the  state  representing  the  strain. 
Hcuristicall>,  then,  one  can  think  of  the  zero  integral  mean  condition 
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ux(t,x)dx  =  o  as  representing  the  boundary  conditions  given  by  (1.2) 


U„ 

Next,  let  X  =  L2(0,1)  with  norm  ||f||^=  j  |f|2dx, 

Y  =  L2(0,1)  with  norm  ||f|£  =  j1  |f|2dx, 
and  W  =  L2(-r,0;Y)  with  norm  ||w|j^,=  J  j|w(s) j|^ds. 


Let  G(s)  and  K  denote  the  2x2  matrices  defined  by 


G(s)  = 


0  0 
o  pg(s) 


and  K 


1  0 
0 


respectively.  Define  the  operator  AQ  :  D(A0)  C  X  x  Y  — ►  X  x  Y 
domain 


D(a0,  =  Hj  (0,1) 


[h1  (0,1)  n  L®(0,1)] 


X 


■DR 


To  obtain  a  wcll-poscd  model,  wc  let  z(t)  denote  the  "state 


sufficient  to  consider  the  norm 


s 


on  Z  . 


The  operator  A  is  not 


dissipative,  and  it  will  be  useful  to  introduce  a  dissipative  operator  Aj  via  a 
similarity  transformation.  To  this  end,  define  the  operator  L  €  E(Z,Z)  by 


(2.3) 

L 

[i] 

LwJ 

.0-W. 

Note  that  L  =  L'1  and  L  is  bijective.  Next  define  the  operator  Aj  on 
the  domain 


D  (Aj)  =  {z  €  Z  :  Lz  €  D (A)} 


.w. 


€  Z 


<P  e  H*(0,1)  ,  w  e  H1(-r,0;Y) 

:  [p1^  ”  p  |°  g(s)  w(s)ds  j  €  H1(0,1) 


w(0)  =  0 


by  Aj  *  L‘!AL  . 


In  particular,  it  follows  that 
[[*♦  “of  g(s)w(s)ds]  ] 
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Well-posedncss  is  provided  by  the  following  result. 


Theorem  2.1.  The  operator  Aj  generates  a  CQ-semigroup. 


Proof.  The  proof  proceeds  as  follows.  We  first  verify  that  Aj  is  dissipative, 

and  then  show  that  RfXI-Aj)  =  Z  for  some  X  >  0  .  We  then  show  that  Aj 

is  a  closed  operator.  It  follows  that  (Xl-Aj)'1  exists  and  is  a  closed,  bounded 

and  densely  defined  operator.  Thus,  DfXI-Aj)'1  =  Z  ,  and  hence 

K(XI-Aj)  =  Z  .  Therefore,  DfAj)  will  be  dense  in  Z  (see  Theorem  1.4.6  in 

[19]),  and  the  Lumer-Phillips  Theorem  will  imply  that  At  generates  a 

|>1 

OsemierouD.  In  order  to  complete  the  details.  let  I  d>  €  D(A,)  .  Then 


C0-semigroup.  In  order  to  complete  the  details,  let  I  ip  €  I^Aj)  .  Then 

Lw  J 

/  I"9!  M  \ 

2  (Aj U  ,  U  )  = 

'  LwJ  LwJ  ' g 

r1  r  f°  l 1  r1 

2  ti/j  -  g(s)w(s)ds  <pdx  +  2  £  (p'tp 

O  L  -r  1  0 

-  2  |  g(s)  |  J^<p '  +  —  j  w  (s)  dx  ds 

f°  f1  dw 

=  -  2  g(s)  —  v.  (s) dx  ds 

-r  J0  dS 

=  ~j  [g(s)  |  w2(s)  dx  j  ds  +  |  w2(s)dxds 

«  f  j-g(s)  [  w2(s)  ds 

f°  f1 

«  H  g0(s)  w2(s)  dx  ds  S  0  . 

•  n 


'-r  ■'0 


Hence,  is  dissipative  in  Z  .  To  show  that  RCXI-Aj)  =  Z  ,  it  is 

sufficient  to  show  that  (XI-Aj)D  =  Z  ,  where  D  is  given  by 


0  =  ^u,  x<p  -  j~0'  +  g(s)w'(s)ds  ^  + 

(v.  ^  ~  ¥')y  + 

(*•  -%)„■ 


D  =  ■  0  e  D (Aj)  :  0  €  H^O.l),  w  €  H^-r.OiH^O.l) )  C  D(Aj)  . 


Suppose  v  is  orthogonal  to  (XI-A)D  .  Then 


T  I 

for  all  0  €  D  .  For  each  <p  €  H2(0,1)  n  Hj(0,l)  ,  let  0  =  y  ip1  and 

1  W  r 

w(s)  =  y  (hence  w(s)  is  a  "constant"  function  in  W  ).  Therefore  using 
€  =  a  +  J°  g(s)ds  j  , 

0  =  (u’  *<P  ~  j5T  “  lp  [y  '?"]  f  g(s)ds  +  j5  f  8(s)  [y«pm]  ds  ^ 

=  <^u,  X<p  -  for  all  <p  €  H2(0,1)  O  H*(0,l)  . 

The  image  of  H2(0,1)  n  Hq(0,1)  under  the  operator  ^X2I  -  |j-2j  is  dense  in 

L2(0,1)  which  implies  that  u  =  0  .  If  u  =  0  and  0  =  0  ,  then 

0  =  (z,  Xw - )  for  all  w  e  H1(-r,0;H1(0,l)) 

\  ds  /w 


-11- 


and  z  =  0  by  similar  reasoning.  Next,  choose  <f  =  0  and  w  =  0  .  It 


follows  that 


0  =  <v,  X0>y  for  all  0  €  H^O.l)  n  Y  . 


This  implies  that  v  =  constant,  and  hence  v  =  0  since  v  e  Y  .  Thus,  we 
conclude  that  R(XI-Aj)  =  Z  .  It  remains  to  be  shown  that  Aj  is  closed. 


Assume  that  the  sequence 


CO 


K 


satisfies 


zn  €  J)(Aj)  , 


z  =  0  in  Z  , 


Aizn  y  =  j  ln  Z 


We  must  show  that  z  €  DfAj)  and  AjZ  =  y  .  We  deduce  from  (2.5)  and  (2.6) 


that  <pn  — ‘if  and 


<P '  — -  8  • 


Since  the  differential  operator  B  ,  defined  by  D(B)  =  Hq(0,1)  and  Bf  =  f1  , 
is  closed  on  L,(0,1)  ,  we  conclude  that 


<f  e  Hq(0,  1) 


m\ 


WWW 
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By  a  similar  argument,  we  see  that 


(2.10) 


g(s)  w(s)ds 


€  H1(0,1) 


and 


(2.10) 


'e  ,  1 

r° 

P  ^  P 

g(s)  w(s)  ds 

« 

-r 

Next,  (2.6)  and  (2.7)  imply  that 
dwn 

— d  —  7-0  =  7-  <p'  m  w  . 
ds 

Also,  (2.5)  implies  that  wn(s)  — *  w(s)  in  W  .  Again,  using  an 
similar  to  that  used  above,  we  conclude  that 


(2.12)  w  €  Hx(-r,0;Y)  , 

(2.13)  w(0)  =  0  , 


and 

dw 

(2.14)  7  =  <p’  +  —  • 

ds 

Combining  (2.8),  (2.10),  (2.12)  and  (2.13),  it  follows  that  z  €  IXAj)  . 
(2.11)  and  (2.14)  imply  that  AjZ  =  y  ,  and  the  result  follows. 


argument 


Also,  (2.9), 

■ 


Observe  that  since  A  =  L‘1A1L  ,  it  follows  from  Theorem  2.1  that  A 


generates  a  C0-semigroup.  We  state  this  in  the  following. 


Corollary  2.1.  The  operator  A  generates  a  CQ-semigroup. 


We  conclude  from  Corollary  2.1  that  the  system  (2.1)-(2.2)  is  well  posed. 


We  are  interested  in  constructing  an  approximation  scheme  for  an  optimal 
control  problem  associated  with  the  viscoelastic  model.  For  S  >  0  and  a 
partition  of  [0,1]  given  by  0  <  Xj  <  1  ,  i  =  1,2, -  - -,P  ,  the  operator 


S®  •*  Lj(0,l) 


is  defined  by 


fXj+6 

sfco  =  i  L  foo* 


The  problem  is  to  minimize  the  cost  functional 


ap  p 

I  |s?  ^(t,x)|2  +  I  |sf  &  (t,x)|2  +  Rju(t)|2  dt 


subject  to  dynamics  governed  by 


~2y(t,x)  =  y(t,x)  +  I  g(s)f-y(t  +  s,x)ds  +  b(x)  u(t), 

t  ox  ox  i  ox 

-r 


where  b(x)  €  L2(0,1)  is  a  given  function.  Initial  data  and  boundary  conditions 
are  as  in  (1.2)-(1.3).  Loosely  speaking,  the  cost  functional  corresponds  to 
observations  of  the  average  velocity  and  strain  in  a  neighborhood  of  each  Xj  .  We 
point  out  that  we  consider  observations  of  average  velocity  and  strain  rather 
than  point  observations  of  velocity  and  strain  because  this  allows  us  to 
formulate  an  LQR  problem  with  bounded  (rather  than  unbounded)  operators  in 
the  cost  functional  [i.e.  point  evaluation  is  not  a  bounded  linear  functional  on 


-  a,  ■  ‘M  fcjj  ‘  >.■  u  ‘W> i^Wi  »  Jj,  t 


Based  on  the  results  of  the  previous  section,  we  can  treat  this  problem 
equivalently  as  an  abstract  regulator  problem  on  the  state  space  Z  : 


minimize 


J8(u)  =  |  jjcsz(t)|2  +  R  |u(t)  |2 j  dt 


subject  to  dynamics  governed  by 


z(t)  =  Az(t)  +  Bu(t) 


z(0)  =  z0  . 


Here  Z  and  A  are  as  defined  in  the  previous  section.  The  operator 
B  :  1R  — ♦  Z  is  defined  by 


b(x) 

Bu  =  0  u(t). 

0 


The  operator  C8  :  Z  — *  IR2p  is  defined  by 


C8  J 


%  .V  *■  -%  /«  -N 


It  can  be  shown  that  (see  [12])  there  exists  a  unique  optimal  control  for  (3.2)-(3.3) 
which  is  given  in  feedback  form  by 

u(t)  =  -Kz(t)  . 

Further,  K  =  R'1B*n  ,  where  n  is  the  unique  non-negative  self-adjoint 
solution  of  the  algebraic  Riccati  equation 

(ARE)  A *n  +  nA  -  rBR_1B*77  +  (C8)*C8  =  0  . 

Since  K  is  a  bounded  linear  functional  on  Z  ,  we  can  write 
u(t)  =  -Kz(t) 

— 

where  kx  €  X  ,  k2  €  Y  ,  and  k3  €  W  .  We  are  interested  in  an  approximation 
scheme  which,  in  addition  to  accurately  simulating  the  dynamics  of  (3.1),  also  gives 
a  reasonable  approximation  of  the  functions  kj  ,  k2  ,  and  k3  . 

We  shall  be  interested  in  approximation  schemes  consisting  of  sequences  of 
finite  dimensional  spaces  ZN  C  Z  and  operators  AN  :  Z  — ►  ZN  ,  each  AN  the 
infinitesimal  generator  of  a  C0-semigroup  TN(t)  .  If  PN  :  Z  — ►  ZN  is  the 
orthogonal  projection  defined  by  ZN  ,  then  for  each  N  we  have  the  following 
finite  dimensional  regulator  problem: 

(3.2)n  minimize  J8(u)  =  j  jjc8zN(t)|2  +  | Ru(t)  | 2  | dt 


subject  to 


ZN(0)  =  PNZn  . 


Here  BN  =  PNB  . 

This  problem  is  finite  dimensional  and  can  be  solved  numerically.  Naturally, 
the  nature  of  the  approximation  scheme  will  determine  how  well  the  solution  of 
this  problem  "approximates"  the  solution  of  the  original  problem  (in  particular, 
whether  the  approximating  feedback  gains  converge).  Gibson  [13]  (see  also 
[4])  studied  this  problem  and  showed  that,  in  addition  to  satisfying  PNz  — ♦  z 
for  all  z  €  Z  ,  an  approximation  scheme  [ZN,AN]  should  have  the  property 


TN(t)z  —  T(t)z  ,  for  all  z 


TN*(t)z  — *  T*(t)z  ,  for  all  z  , 


with  uniform  convergence  on  bounded  t-intervals.  Here,  T(t)  denotes  the 
semigroup  generated  by  A  and  T*(t)  (=  T(t)* )  is  the  semigroup  generated 

by  A  We  note  that  an  approximation  scheme  (ZN,AN)  which  satisfies  only 

(3.4)  is  useful  for  "simulation"  purposes;  i.e.,  approximation  of  the  open  loop 
problem.  However,  for  the  closed  loop  problem,  (3.5)  is  important  in  order  to 
guarantee  strong  convergence  of  the  approximating  feedback  gains.  For  the 
approximation  scheme  which  we  shall  develop,  convergence  results  (i.e.  (3.4)  and 

(3.5) )  are  given  in  the  context  of  following  version  of  the  Trotter-Kalo 
semigroup  approximation  theorem  (see  [19]). 


>./•/, -\ /,/. v  v.  .‘,v. /.  «*. 


Theorem  3.1.  Let  A  €  G(M,0)  ie  //re  infinitesimal  generator  of  a  CQ-semigroup 
T(t)  in  a  Hilbert  space  Z  ,  and  suppose  there  is  a  sequence  of  linear  operators  AN 


each  of  which  generates  a  CQ-semigroup  on  Z  .  If 
HI)  An  e  G(M,0)  for  N  =  1,2,- ••  . 


H2)  Anz  — »  Az  for  z  €  D  ,  D  a  dense  subset  of  Z  , 


and  H3)  there  exists  with  Re  XQ  >  0  such  that 


(A  -  X0)D  is  a  dense  subset  of  z. 


then  TN(t)z  — *  T(t)z  for  all  z  e  Z  ,  t  >  0  ,  and  the  convergence  is  uniform  in 


compact  t-intervals. 


With  this  preliminary  discussion  in  mind  (especially  the  importance  of 


(3.4)  and  (3.5)),  we  proceed  to  develop  an  approximation  scheme  for  our 


problem.  Recall  that  the  state  space  Z  is  defined  by 


Z  =  XxYxW  =  L,(0, 1 )  x  L°(0,1)  x  L2{-r,0;L°(0,l)) 


with  norm 


<p  2  »i  f°  f1 

<l>  =  ||<p|l£  +  II^Hy  +  ||w||w=  ((f2+4'2)  +  w2(s) dx  ds 

\i/  •  c\  n 


The  state  operator  A  is  defined  on  the  domain 


<f  €  Hq(0,1  )  ,  w  €  HH-r.OjY) 


3(A)  -  U  €  Z  : 


+  ^  g(s)  w(s) dsj  €  HH0.1)  ,  w 


(0)  =  <p 


'j-  "j-  's-  s 


9 

A  ip 


^  [p*  +  p  J  g(s)w(s)dsj 


ao{k[I]  +  f  G(s)  [w°(s)]  ds 


We  will  also  need  to  consider  the  adjoint  of  A  .  It  is  straightforward  to 
verify  that  A*  is  defined  on  the  domain 


(3.10)  D(A  )  =  U  €  Z 


if  €  Hj(0,l)  ,  iP  e  HHO,!) 


w  €  H'l-rAY)  w(-r)  =  0 


(3.11) 


A  <P 


=  “  o'?'  +  w(°) 


-w  -  -  g(s)<p' 


Ao  "p  0  J  "  0 

lo  -iJl?J  Lw(o) 


dw  1 

-  -  -  p  (!s),' 


Equation  (3.9)  suggests  a  two-stage  approximation  for  A  .  In  particular,  one 


must  approximate  the  operator  A0  by  discretizing  the  spatial  variable  x  . 

d 

Further,  one  must  approximate  —  by  discretizing  the  delay  variable.  Roughly 

ds 

speaking,  we  will  use  a  finite  element  approximation  of  A0  ,  and  the  AVE 

approximation  scheme  (sec  [2])  for  the  delay  variable. 

To  proceed,  divide  [0,1]  into  N  equal  subintervals  and  let  hj^(x)  , 

i  =  0,1,- ■ -,N  ,  denote  the  standard  "hat"  functions  which  are  continuous, 

piecewise  linear,  and  satisfy  h^(j/N)  =  6^  ,  for  j  =  0,1,  •  • -,N  .  For 

convenience  in  the  ensuing  discussion,  we  shall  use  the  shorthand  notation 

h  =  hN(x)  and  h1  =  —  h^(x)  .  Next,  we  define  the  subspace  XN  of  X  bv 
1  1  1  dx  * 

XN  =  span|h;|i  =  1 ,2,  •  •  •  ,N- 1 } 

and  the  subspace  YN  of  Y  by 

Yn  =  span  {V1  |i  =  1,2,.  ,N-1 )  . 

Next,  define  the  bilinear  form  a(-,  )  on  X  x  Y  by 


We  define  the  operator  Aq  (the  finite  element  approximation  of  A0)  as  the 
restriction  of  a(  - ,  - )  to  XN  *  YN  : 

(3.13)  <A*j  u,v  >  =  a(u,v) 

u  X*Y 

for  u,v  €  XN  x  YN  . 

To  complete  the  "first  stage"  of  the  approximation,  let 


WN  =  L2(-r,0;YN)  C  W 

(3.14) 

ZN  =  XN  X  \'N  X  wN  , 


and  define  AN  on  the  domain 


denote  the  orthogonal  projection  onto  ZN  ,  and  note  that 


where  :  X  — ♦  XN  ,  Py  ;  Y  — *  YN  ,  P^.  :  W  — »  WN  arc  the  respective 

orthogonal  projections. 

Observe  that  for  each  N  ,  AN  is  the  infinitesimal  generator  of  the  C0- 
semigroup  on  ZN  corresponding  to  the  ordinary  delay-differential  equation 


(3.18) 


Thus,  for  the  "second  stage"  of  approximation,  we  can  use  any  of  several 


schemes  in  the  literature  for  this  type  of  problem  (sec  [2],  [3],  [17]).  Since  wc 
are  interested  in  optimal  control,  in  this  paper  we  will  use  the  AVE  scheme  (sec 
[2]),  which  is  known  to  have  desirable  properties  (recall  (3.4)  and  (3.5))  for  such 
problems  (see  [13]). 

Subdivide  [-r,0]  into  M  equal  subintervals  [t^,t^]  ,  j  =  1,2,  ■  ,M  , 
where  t^1  =  -  .  Let  X^V)  denote  the  characteristic  function  of  [tM.tMj  , 

[tj^t^j]  ,  and  define  the  spaces  WN,M  C  WN  and  ZN,M  C  ZN  by 


(3.19) 


W 


•  N.M 


w  ‘ 


w- 


M 

r,0;Y)  |  wN  =  E  v;Nx^(-), 
j=i 


vN 

j 


e  Yr 


and 


(3.20) 


ZNM  =  XN  x  YN  X  WN'M 


respectively.  Next  define  the  operators  QNM  :  zN  — *  ZN,M  by 


where 


-f 


•J!. 


(3.22) 


(wN)M  =  0N,  (wN)M 


M 


wN(s)  ds 


UNUMCUHjCUh 


^  ^  v  \. 


-V  »w  V  A.  V  ».  S.  ■-  V' 


(3-23) 

w  M 
GM  =  — 

J  r 

f  J  1  G(s)ds  , 

tM 

j 

and  gM  = 

M 

r 

for  j  = 

1,2,  ■  ■  ,M  . 

Combining  the 

two  stages 

of 

operator 

Anm  :  Z 

_  ZN,M  by 

(3.24) 

AN’M  =  Q 

N,M  pN 

rz  • 

g(s)ds  , 


Remark  3.1.  In  view  of  hypothesis  (H2)  of  Theorem  3.1,  the  approximation 
scheme  should  satisfy 

||AN,Mz  _  Az||  _  o 

for  all  z  in  some  dense  subset  of  Z  .  From  the  triangle  inequality, 

||An,mz  -  Az ||  S  ||Qn,mPnz  -  AnPnz||  +  ||ANPNz  -  Az||  =  Sj  +  S2  . 

Standard  spline  estimates  imply  that  S2  — ►  0  as  N  — *  00  for  each  z  ,  but 

convergence  of  the  AVE  scheme  implies  only  that  Sj  — *  0  as  M  — ►  00  for 

each  fixed  N  and  z  .  In  particular,  the  rate  of  convergence  in  M  of 

qN.m  tQ  an  js  bounded  by  ||Aq  j|  .  Although  A0  is  an  unbounded 
operator,  it  is  straightforward  to  show  that  ||Aq||  =  O(N)  (sec  Lemma  3.2). 
This  estimate  provides  the  key  to  the  choice  of  the  index  M  as  a  function  of 
the  index  N  so  that  the  convergence  in  M  dominates  the  unbounded 

behavior  of  ||Ao||  (see  3.43). 

We  now  state  some  useful  properties  of  the  operator  Ag  and  the 
projections  given  in  (3.17). 


I**  *'■.  /.  /.  /.  r,  .*  /,  -"...V.V,  -V-'\  O 


Lemma  3.1. 


i) 

If  f  €  X  ,  then 

(3.25) 

iip!!f  -  rnx  - 0 

as 

N 

— *■  CD 

ii) 

If  f  €  H^O,!)  C 

X  . 

then 

(3.26) 

l!px<r'>-  f,llx-  0 

as 

N  —  » 

and 

(3.27) 

||<Pxf>'  -  f  iu  - 0 

as 

N  —  oo 

in) 

If  f  €  Y  ,  then 

(3.28) 

llp?r  -  fllY  ~ 0 

as 

N 

- ►  CO 

iv) 

(3.29) 

PV(S)  =  PyW(s) 

in 

the 

L2  sense. 

v) 

If  w  €  W  .  then 

(3.30) 

llP>  “  Wllw- 

8 

I 

Z 

O 

vi) 

//  [ l]  €  D(A0)  . 

then 
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vii)  If  we  H1(-r,0;Y)  C  W  ,  then  P^w  e  H^-r.OiY*)  C  WN  .  and 

(M2)  p*  K )  °  J7  p*,w(s) 


Proof.  For  i)  and  ii),  one  can  use  virtually  the  same  argument  used  by  Banks 
and  Kappel  in  ([3],  Theorem  4.1),  with  the  slight  modification  that  our  splines 
satisfy  zero  boundary  conditions.  For  iii),  let  f  €  Y  =  Lj(0,l)  and  define 

F(x)  =  J*f(T)dT  .  Then  since  Py  is  the  orthogonal  projection  of  Y  onto 
Yn  ,  it  follows  that 


min 

u€YN 


f 


Y 


5  ll<pxF>'  -  F'llv-  0 

by  (3.27).  To  prove  iv),  we  note  that  it  follows  from  the  definition  of  Py  as 
an  orthogonal  projection  that 

0  =  [  [Py[w(s)]  -  [w(s)])-cN(x)dx  for  almost  all  s  €  (- r,0)  , 
for  all  cN(x)  €  Yn  .  Hence,  it  follows  that 

0  =  J°  |  (Py[w(s)]  -  [w(s)]j  bN(s)dx  ds 

J  -r  J  a 

for  all  bN(s)  €  WN  .  Now  iv)  follows  since  P^,  is  the  orthogonal  projection 

of  W  onto  WN  .  Next,  v)  follows  from  (3.28),  (3.29)  and  the  dominated 

convergence  theorem.  To  see  that  (3.31)  holds,  assume  that  [$]  e  D(A0)  and 


that 


Finally,  to  verify  that  vii)  is  true,  let  w  e  H^-r.OiY)  C  W  .  Then 

w(s)  =  w(-r)  +  Jsr  w(T)dr  .  Using  (3.29)  we  have,  for  all  c(s)  €  W  , 

0  =  ^P>(s)  -  P^jw(-r)  +  |  w(T)dT]  -  c(s)  ) 

=  ^w(s)  -  ^w(-r)  +  P^  J*  w(T)dT  ,  c(s)  ^ 

Since  Py  is  bounded,  it  can  be  moved  inside  the  integral  (see  [8],  p.  91). 


Hence, 


=  ^P^w(s)  -  [p^w(-r)  +  J  P^w(T)dTj  ,  c(s)  ^ 

=  (p>(s)  -  [p^w(-r)  +  |  pN,w(T)drj  ,  c(s)  ^  . 

Thus,  P^w(s)  =  P^,w(-r)  +  J“rP^w(T)dT  ,  which  implies  that 

Pww(s)  6  HH-rAY)  and  ~  P^w(s)  =  ?w  ^  • 


Thus,  P^w(s)  =  P^,w(-r)  +  J®rP^w(T)dT  ,  which  implies  that 


P vyW(s)  €  H!(-r,0;Y)  and  P”w(s)  =  pw  ^7 


The  next  lemma  provides  the  growth  rate  for  IIAqII  . 


Lemma  3.2.  There  exists  a  constant  K  ,  independent  of  N  and  uN  ,  so  that 


<3-38>  !IaouNI!x*y  4  kn2I!uNIIxxy 


for  all  uN  e  XN  x  Yn  . 


where  Uj  =  E  a jhj  and  u2  =  E  b^'  .  Let 


fu  1 

Proof.  Let  uN  =  1  ,  where  ut  =  E  a  jhj  and  u2  =  E  bjhj1  .  I 

vUj  J  i  =1  j  =1 

fu  1  fv  1  N’X  N_1 

Aq  1  =  1  ,  where  Vj  =  E  orh;  and  v2  =  E  Bjh’  .  By  a  direct 

lu2J  Lv2J  i=l  i=i 

calculation  using  the  definition  of  we  have 


N-l  N-l 

(3.39)  E  otj<hi,hj>x  =  -  E  b^h/,^' >x 

i=l  i=l 


N-l  N-l 

(3.40)  E  B; <h ,hj'  >Y  =  E  a|<h.',hj'>Y  , 

i=l  i=l 

for  each  j  =  I,2,-  -,N-1  .  Hence,  using  (3.40)  twice  we  have  that 


(3.41) 


2  /  N-l  N-l  v 

(  z  6,h,\  z  8,1.;) 

Y  '  i=l  ;=1  'Y 

/  N-l  N-l  v 

=  (  e  ®iV.  ^  ^ihi7 

'  1=1  i=l  'Y 


N-l 

E  aih' 

i=l 


<  K-N2-  E1aihj  =  KN2|JUl , 

i=i  y 


where  the  last  inequality  follows  by  applying  the  Schmidt  inequality  ((20])  over 
each  interval.  Also,  by  using  (3.39)  and  the  Cauchy-Schwartz  inequality,  we 
have  that 


VA  -.  AA 


J  N-l  N-X  y 

(  Z  aihi’  1  «iO 
'  i=l  i=l  'X 

/  N-l  N-l  y 

-  (  £  b,h’,  z  «,!,,) 

i=i  i=i  'x 


s  wx  r;nx « ',k  Niiuiiix'iiviiix  • 


where  the  last  inequality  follows  by  again  applying  the  Schmidt  inequality  over 
each  interval.  Hence, 


(3.42)  v.  S  KN2  u. 

II  xilx  l  x 


The  result  follows  from  (3.41)  and  (3.42). 


We  are  now  in  position  to  prove  the  necessary  convergence  results  for  the 
approximation  scheme  defined  above  in  (3.20)  and  (3.24).  As  discussed  above  in 
Remark  3.1,  it  will  be  necessary  to  impose  the  following  condition  on  the 
indices  N  and  M  . 


Definition.  A  sequence  (N,Mn)  satisfies  condition  Cy  if  for  all 


N  =  1,2,-  •  •  , 


(3.43) 


*°  Mn  mn  mn 

£  gj  Xj  N  -  g(s) 
-r  3=1 


ds  <  ^7’  ^  >  1  - 


where  the  constant  K  is  independent  of  N  . 


.N  .  -  .V, 


We  note  that  due  to  the  convergence  properties  of  the  AVE  scheme,  a 
sequence  satisfying  condition  Cy  exists  for  any  g(s)  6  L2(-r,0)  .  In  particular, 
for  g  €  C1(-r,0)  ,  the  condition  is  satisfied  by  choosing  MN  =  N?  ,  7  >  1  . 

Given  a  sequence  (N,Mn)  ,  let 

(3.44)  =  AN’Mn  , 

and  denote  by  TN(t)  the  C0-semigroup  generated  by  /1N  on  Z  . 


Lemma  3.3.  (Stability)  Let  (N,Mn)  be  a  sequence  satisfying  condition  Cy 

Then  there  exists  constants  such  that  for  all  N  £  1  , 


TN(t)  S  Me0t. 
z 


Proof.  For  each  M  ,  define  the  norm 


on  Z  by 


<p  2  ft  pO  M  fb 

0  =  (p<p2  +  € 02)  -  [ig^X*4]  w2(s) 

W  J0  J.r  j=l  Ja 


dx  ds  . 


These  norms  are  uniformly  equivalent  to  the  original  unweighted  norm  on  Z 
It  is  sufficient  to  show  that  L‘1AN,ML  is  dissipative  in  the  II  -II  M  norrr 


To  show  this,  let  z  =  0  €  Z  .  Then 


<L'1An-mLz,z>m  =  <QN’MPNLz,L*z> 


where  L  is  the  adjoint  of  L  with  respect  to  the  gM  norm  and  is  given  by 


r<pT  j  »o  m 

L*  0  =  0  -  j  (  E  gf* xf )  w 

w  '  -r  j=l 


Therefore, 


♦  “  f'  [(»■")»  -  (w»)M]gM(wN)Mdx  . 

j=X  ^0 


The  first  two  lines  cancel  exactly,  leaving  only  the  term 
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M  ,.1 

E  P)M  -  (wN)M]gM(wN)Mdx  £  E, 
j=l  0  '  J  J 

Using  Cauchy-Schwarz  and  the  inequality  2ab  <  a2  +  b2  ,  we  get 

E,  «  “j  g?}’  {[<«">»  ]’  -  [(w^fjdx 

-  J  g“  {  [(wN)M]!ds  -  ;  g|j  |  [lwN')“],d.  ♦  Z  (g“ ,  -  gf)  [  [(wN)M]!dx 

J0  Z  J0  j=l  z  Jo 

5  \  8^1  [(wTfdx  =  \  8^l  (^N)2  ■ 


1  • 


The  last  inequality  follows  since  g  is  increasing,  and  hence  gJ*j  S 


Finally,  because  g  is  bounded,  the  gy1  are  bounded  iniformly  in  M  .  We 
conclude  that 


<L_1AN,MLz,z>g  «  K<z,z>g 


uniformly  in  N  and  M  . 

Next,  define  the  set  D  C  0(A)  by 


D  = 

'  <P' 

€  0(A)  :  w  e  C*(-r,0;Y) 

.W. 

4 

Lemma  3.4.  (Consistency)  Let  (N,Mn)  be  a  sequence  satisfying  condition  Cy 

n,mn 

and  such  that  MN  — ►  00  as  N  — *  00  .  Then  A  z  — *  Az  as  N  — ♦  00  for 


all  z  €  D  . 


Proof.  Let  z  = 

have 


e  D.  Then,  setting  M  =  MN  for  notatirmal  ease,  we 


I  n,m,, 

A  Nz  -  Az  = 


an  • 
Ao 

K 

r  n 
<P 

i.  • 

M 


+  I  ~G; 


M 


j=l 


M  J 


r 

o 

' 

r<£>3  f°  fn  1 

N 

(WN)M 

V  J 

“  Ao' 

K 

y*Ki 

ds 

XxY 


r 

M 

I 

m  r 

rpjwiM  -  iwflxf  -  - 

i. 

j=l 

T  L 

L  w  Jj.j  L  w  -lj  J  J  ds 

ds 


=  Sj  +  S2  . 


Estimating  the  term  St  ,  we  have 

.0  M 


Si  < 


[  A"  E  [c“x«  -  G(s)]  [  l 

’  -r  j  =  l  v 


0  1 

(s) 


ds 


XxY 


K  (£]  +  fG(s>  [„°(s)] ds}  -  A»{K  (*]  *  frG(s) 


ds 


XXY 


"  F1  +  F2  • 


It  follows  from  (3.31)  of  Lemma  3.1  that  F2  — *  0  .  For  the  term  Fj  ,  apply 
first  the  Cauchy-Schwarz  inequality  and  then  Lemma  3.2  to  get 

Fi  «  f  I  |gjMxf  -  g(s)|2ds.N2- ||wN||* 

J-rj=l 

Since  N,Mn  satisfy  condition  Cy  ,  it  follows  (see  3.43)  that  Fj  — *  0  . 


Finally,  considering  the  term  S2  ,  we  have  the  estimate 


w  *  T  K  •]“  -  -  d7  K  »)  ds 

-I*  j—  1  y 


d  M  dw 

—  PN  w - ds 

ds  w  ds 


=  Ej  +  E2  . 


It  follows  from  (3.30)  and  (3.32)  of  Lemma  3.1  that  E2  — ♦  0  .  For  the  term 
Ej  ,  wc  follow  an  argument  analogous  to  that  used  by  Banks  and  Burns  in  the 
proof  of  Corollary  3.1  in  [2],  This  leads  to  the  estimate 


where 


E  5  r.f  sup  i  Mj2  +  _L  fl  M  +  IS]' 

1  llSjSM  2  ■>  J  M  12  1  2  J 


r 

rf  •  sup  ||(ip>](9)  -  [^Pj](T)|  :6.t  e 

>.  Y 


I  lid w  dw  r  \a 

‘  sup  :e-T 


K  =  sup 


(6)  :0€[-r,O] 


Since  N  — ♦  00  and  M  — *  00  ,  it  follows  from  the  uniform  continuity  of  — 

ds 

on  [-r,0]  that  Ej  — k  0  . 


Lemma  3.5.  For  the  set  D  defined  above,  there  exists  a  real  number  y0  such 
that  (A-XI)D  is  dense  in  Z  for  all  Re  X  >  7n  . 


-35- 


Proof.  Since  A  generates  a  C0-scmigroup,  there  exists  J0  such  that 


Re  \  >  yQ  implies  X  €  p(A)  .  Fix  X  with  Re  X  >  yQ  .  For 
the  equation  (A-X) 

dense  subset  S  C  Z  by  S  =  XxY  x  C(-r,0;Y)  .  If 
satisfies 


If1 

'a' 

f<p  j 

= 

B 

has  a  unique  solution 

$ 

Lw  J 

.z  J 

lw  J 

a 

B 

iz  J 


6  Z  . 


B 

iz  J 


€  IXA)  .  Define  the 
€  S  .  then 


(3.45) 


a  1  f° 

p  <l>  +  p  I  g(s)w(s)ds 


—  X<p  =  a 


(3.46) 


<p  1  -  >  i =  B 


d  w 

(3.47)  —  -  Xw  =  z  . 

ds 

Since  w  e  H1(-r,0;Y)  and  z  €  C(-r,0;Y)  ,  it  follows  from  (3.47)  that 
w  e  C^-dOiY)  .  Hence,  S  C  (A-X)D  and  the  result  follows. 

■ 


The  following  is  an  immediate  consequence  of  the  above  results. 

Theorem  3.2.  Let  N,Mn  be  a  sequence  satisfying  condition  .  and  such  that 

Mn  — *  00  as  N  — ♦  ®  .  If  the  sequence  of  operators  AN  :  Z  — *  ZN,M  is 

defined  as  above  (sec  (3.44)),  then  the  hypotheses  of  the  Trotter-Kato  theorem  hold; 
that  is.  TN(t)  —  T(t)  strongly  on  Z  . 

■ 

Next,  we  verify  a  similar  result  for  the  adjoint  A  (sec  (3.10)  and 
(3.11)).  Since  |TN,M(t)||  =  ||TN,M(t)*||  ,  we  need  only  verify  hypotheses  H2)  and 
H3)  of  the  Trotter-Kato  theorem.  Recall  that 


r  ’j*  V  V  V 


AN'M  =  QN'MPN  =  QN,MnMpN  _  pNnMQN,MnMpN 

where  nM  =  nN,M  :  ZN  — *  ZN,M  ;s  the  orthogonal  projection. 

It  is  useful  to  distinguish  between  QN,M  (defined  on  all  of  ZN)  and  its 
restriction  to  ZN'M  .  Let  RN'M  =  QN'M  ^  N  M  :  ZNM  —  ZN'M  .  Hence, 

AN,M  _  pNpjMj^N.Mj^MpN 


aN,M*  _  pNnM(RN>M)‘nMPN 

where  the  adjoint  is  computed  in  the  unweighted  inner  project  (so  that  the 
projections  PN  ,  nM  are  orthogonal).  To  determine  the  action  of  AN,M’  ,  we 
compute  matrix  representations  relative  to  the  following  basis.  For  i  =  I ,  •  •  •  ,N- 1 
and  j  *  2,  •  ■  -  ,M+1  ,  let 


ei  =  o 


e(N-l)  +  i 


:  =  h  > 


'j(N-l)  +  i 


The  matrix  representation  of  (RN,M)*  is  given  by 


M'i(RN'M)TM  , 


where  M 


is  the  (M+2)(N-1)  *  (M+2)(N-1)  matrix  whose  i-j 


entry  is  <e;,ej>z  >  and  rn- M  is  the  matrix  representation  of  RNM  .  Here,  the 
superscript  "T"  denotes  the  matrix  transpose.  RN* M  is  calculated  according  to 
the  formula 

=  M-1[<RN'Mei,ej>z]T  . 

We  have 


D 


1 


is  also  an  (N-l)x(N-l)  matrix  . 


Next,  define  the  set  D  C  D(A  )  by 


D*  =  U  €  D( A*)  :  w  €  C1  (-r,0;Y) 


If  4>  €  D  ,  then  from  the  above  matrix  representations  it  is  clear  that 


-p  o  <PN 

,  N  P  + 


0  -lJUNJ  l(wN)* 


(An'm)  0  = 

.w  J 

I  M  M  v* 

•^N  r  -  i  "  [(wN)»-(w")“,]x» 

j=i  j=t 


where  (wN)m+,  is  defined  by  (wN)m+i  *  w(-r)  =  0  (compare  with  (3.11)). 


Lemma  3.6.  Let  (N,Mn)  be  a  sequence  satisfying  condition  Cy  .  and  such  that 
Mn  — *  ®  as  N  — *  00  .  Then  (AN,M)  — ♦  A*z  as  N  — *  00  for  all  z  €  D 


A 
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Proof.  The  proof  is  quite  similar  to  that  of  Lemma  3.4. 
then 


AK'M)*z  -  A*z  2  < 


an 

Ao 

a  N 
-p<p 

o 

< 

1 

f  a 
-p<P 

.  -*N  J 

.  -<1>N 

X*Y 


+  ![.*]“  -  w(o) 


f° 

M 

1 

-  p  <p'  g(s) 

i  gjM  xf* 

J  -r 

j=i 

♦f 

M  M 
I  f 

[(wN)M  _ 

(wN)M+1]XM 

*  -r 

j=0 

ds 


ds 
d  w 


ds 


-  Fj  +  F2  +  Fg  +  F4  . 


Fj  — •  0  and  F4  — *  0  as  before.  For  the  term  F2  , 


F2  *  [wn]M  -  wN(0)  2  +  wN(0)  -  w(0) 

1  V 


-  Sj  +  s2  . 


Lemma  3.1  implies  S,  — ►  0  .  Also, 


Si  < 


KA 


2r2 


M2 


where  K  =  sup{|jw(6)  :  8  e  [-r,0]  } 


(use  (3.22)  and  Cor.  3.1  in  [2]).  Hence,  F,  — *  0 


Let 


£  -e-  -6 


Mil  MlMIfe'lMlMlM  'j 


Finally, 


f°  i  i  r  m  i 2 

P  <fN  £  gjM  XjM  -  g(s)  ds 

J-r  L=1  J 

Y 

10  1  i  2 

+  p  [<PN '  —  <f]  g(s)  ds 


“  Ei  +  E2 


Clearly  E2  — »  0  ,  and  it  follows  from  the  Schmidt  inequality  (see  the  proof  of 
Lemma  3.2)  and  condition  Cy  that  Ex  — *  0  .  The  result  follows. 


Lemma  3.7.  For  the  set  D  defined  above,  there  exists  a  real  number  yQ  such 
that  (A  -  XI)D*  is  dense  in  Z  for  all  Re  X  >  yn  . 


Proof.  Use  an  argument  similar  to  that  used  in  the  proof  of  Lemma  3.5. 


Theorem  3.3.  Let  N,Mn  satisfy  condition  Cy  .  Then  TN,M(t)*z  — *  T(t)*z  for 


all  z  €  Z  uniformly  on  compact  X-intervals. 


Therefore,  from  Theorem  3.2  and  Theorem  3.3,  we  sec  that  the 
approximation  scheme  defined  above  satisfies  (3.4)  and  (3.5).  As  discussed 
above,  we  conclude  that  this  scheme  provides  a  reasonable  approximation  for 
simulation  purposes  (the  open-loop  problem)  as  well  as  for  purposes  of 
constructing  feedback  gains  for  the  (closed-loop)  regulator  problem.  In  the  next 
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section,  we  give  the  results  of  some  numerical  experiments  using  this 
approximation  scheme. 


4.  Numerical  Results 


In  this  section  we  present  some  numerical  results  related  to  the 


approximation  scheme  which  we  have  discussed.  Let  gp(s)  denote  the  function 


gp(s)  = 


e5s(-s) 


M 


-e5a[r(s  +  £)  +  (-f)^ 
where  the  constants  f  and  p  are  related  by 

1  +  (e-5r(p/r)^/5 


■  «  s  ?  -r/p 
r/p  S  s  «  0, 


(4.1) 


r  = 


L  ((e-5r/p-l)  +  (5r/p))/25 


(4.2) 


and  1  <  p  <  +  Condition  (4.2)  implies  that 


r° 

gp(s)ds  <  I  and  for  the 


parameters  p  =  a  =  1  used  in  the  numerical  runs  below,  it  can  be  shown  that 


gp  satisfies  Conditions  (1)  and  (2)  in  Section  1  above.  Observe  that  as  p 


g  (s)  -  -c5o(-s)'^  and  gj(s)  is  a  linear  function.  The  selection  of  the  form  of 


the  function  gp(s)  is  made  to  insure  that  Conditions  (1)  and  (2)  hold  and  to 


investigate  the  case  where  gp(s) 


-(-s  r5*  at 


s  =  0.  This  singular  case  is 

important  in  damping  studies.  In  particular,  it  can  be  shown  that  if  g(s)  = 

-e5’(-s)'^,  then  the  (open-loop)  eigenvalues  of  A  are  asymptotic  to  a  quadratic 

curve  and  if  g  e  H1(-r,0),  then  the  (open-loop)  eigenvalues  of  A  are 

asymptotic  to  a  vertical  line  (see  [14]  for  details).  Although  g  (s)  belongs  to 


HV^O),  for  p  =  210  the  function  gp(s)  is  "numerically  singular"  and  the 
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numerical  results  for  this  case  should  be  indicative  of  a  truly  singular  kernel. 
Recall  that  for  g(s)  s  0,  the  undamped  system  is  the  wave  equation  with 
fundamental  frequencies  =  krr,  k  =  1,2, •••  . 

Example  4.1.  For  this  run  we  set  p  =  210  and  constructed  the  approximate 
operator  AN,M.  The  IMSL  routine  ETGRF  was  used  to  compute  the  eigenvalues 
of  An,m  for  various  values  of  M.  Since  we  are  interested  in  the  damping 
properties,  we  display  only  those  eigenvalues  XN’M(k),  k  =  1,2,  ••■N-l, 
corresponding  to  the  first  N-l  fundamental  frequencies. 

Figure  1  illustrates  the  behavior  of  X8,M(k),  for  M  =  4,8,16,32  and  64. 
The  interesting  feature  here  is  that  for  low  values  of  M  the  damping  curve 
predicts  near  viscous  damping  and  as  M  increases  the  damping  curve  becomes 
quadratic  (as  to  be  expected  for  "singular  kernels").  This  figure  supports  the 

remark  made  earlier  that  condition  Cy  is  indeed  necessary.  We  made  several 

other  runs  for  other  values  of  N  with  precisely  the  same  qualitative  results. 

Example  4.2.  We  next  consider  the  optimal  control  problem  ((3.2)N,  (3.3)N).  For 
this  example  we  used  the  kernel  g(s)  =  gp(s)/5  with  p  =  210,  and  selected  the 
observation  points  at  x(  =  .25,  .32,  .50,  and  .67,  i  =  1,2, 3, 4.  As  before,  a  = 
p  =  1  and  b(x)  =  x2,  0  <  x  <  1.  Potters  method  was  employed  to  compute 
the  optimal  feedback  gains.  Shown  in  Figure  2  are  the  open  loop  poles  for 
N  =  8  and  M  =  4,  8,  12,  16,  and  32.  The  solid  dots  are  the  closed  loop  poles 
for  M  =  32. 

Example  4.3.  Using  the  same  data  as  in  the  previous  example,  we  next  examine 

the  convergence  of  the  feedback  gains.  Recall  that  the  feedback  gain  for  the 


rft.»*fcr| 


infinite  dimensional  problem  (3.2).  (3.3)  can  be  characterized  hv  funct:ons 
kj  €  X,  k2  e  Y  and  k3  €  W.  Similarly,  the  gain  for  the  approximating 
problem  (3.2)N,  (3.3)N  can  be  characterized  by  functions  kj*  €  XN,  k^  €  YN, 
and  k3,M  €  WN,M.  In  Figure  3  we  observe  convergence  of  the  functions  k^(x) 
as  N  and  M  increase.  We  used  values  of  M  much  larger  that  those  of 
N  in  order  to  satisfy  condition  Cy.  In  Figure  4  and  Figure  5,  we  observe 
similar  convergence  for  k^  and  k3’M. 


We  conclude  with  the  following  remarks.  Due  to  the  size  of  the 
approximating  system  (for  given  N  and  M,  N  x  (M+2)  is  the  dimension  of 
ZN,M),  and  the  requirements  on  the  discretization  parameters  N  and  M  (i.c. 
condition  Cy),  the  time  and  storage  requirements  for  computation  quickly 
become  unreasonable.  Hence  these  numerical  results  are  preliminary  in  the  sense 
that  we  could  not  "push"  the  scheme  much  beyond  N  =  10.  We  are  currently 
working  on  a  related  unconditionally  convergent  scheme  (i.e.  no  condition  like 
Cy)  with  the  hope  of  reducing  this  computational  burden.  We  are  also 
investigating  the  applicability  of  these  schemes  to  fractional-power  damping 
models  in  beams  (sec  [1],  [20]). 
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